Fermi surface evolution in the antiferromagnetic state for the electron-doped t-t'-t"-J 

model 
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By use of the slave-boson mean-field approach, we have studied the electron-doped t-t'-t"-J model 
in the antiferromagnetic (AF) state. It is found that at low doping the Fermi surface (FS) pockets 
appear around (±7r,0) and (0, ±7r), and upon increasing doping the other ones will form around 
(±-|,±^). The evolution of the FS with doping as well as the calculated spectral weight are 
consistent with the experimental results. 

PACS numbers: 74.25.Jb, 71.10.Fd, 74.72.-h, 74.25.Ha 



Recent angle-resolved photoemission spectroscopy 
(ARPES) measurements^"^ have revealed the doping evo- 
lution of the Fermi surface (FS) in electron-doped cuprate 
Nd2-a:Cea;Cu04. It was found that at low doping a small 
FS pocket appears around (tt, 0), in contrast to the hole- 
doped case where the low-lying states are centered at mo- 
mentum {j, j). Upon increasing doping another pocket 
begins to form around (§, f) and eventually at optimal 
doping X = 0.15 the several FS pieces evolve into a large 
curve around {tt,tt). 

The ARPES data provide clear evidence for electron- 
hole asymmetry in high-temperature superconductors. 
This asymmetry has already been seen in the tempera- 
ture/doping phase diagram, where the antiferromagnetic 
(AF) phase is much more robust and at the same time 
the superconducting (SC) one is much narrower in the 
electron-doped materials. To understand it, the role of 
the long-range hoppings t' and t" was stressed. Based on 
the t-t'{-t")-J model, the single-hole/electron behavior 
and its excitation spectrum have been studied by various 
methodsi^i^i^iSii It was shown that a doped hole enters 
into the lowest energy state at momentum (f , f ), while 
a doped electron goes into that at (7r,0). The excitation 
spectrum for slightly hole/electron-doped t-t'-t"-J model 
has also been investigated by variational Monte Carlo 
technique^ which provides an explanation for the differ- 
ent FS structure in the hole- and electron-doped cuprates 
at low doping. In all these studies, the electron-doped 
Hamiltonian has been mapped onto the hole-doped one 
with extra minus signs for the hopping parameterSf^iS 
thus both electron- and hole-doped systems have been 
treated in the same manner. 

On the other hand, it was claimed recently by Kusko 
et al^ that the use of t-U type models is essential in the 
study of electron-doped cuprates. By mean-field (MF) 
treatment on the t-t'-t"-U model they have found that 



the FS evolution with doping in the AF state is consistent 
with the ARPES results. But in their theory the on-site 
U is treated as a doping-dependent effective parameter 
phenomenologically. 

It is not clear whether a different Hamiltonian should 
be constructed when electron-doped cuprates are stud- 
ied. At least one would ask if the t-t'-t"-J model could 
have reproduced the FS evolution in the electron-doped 
case (even without argument of any doping-dependent 
parameter as done by Kusko et ai). As already stated 
above, for single or slightly electron-doped case the t- 
t'-t"-J model has given good results compared with the 
experiments. Naturally the same model should be further 
studied at finite doping. 

Another concern for the choice of t-J type models is 
the SC pairing symmetry. It is now generally believed 
that the SC gap in hole-doped cuprates has a d-wave 
symmetry. And the theoretical studies were substantially 
based on the t-J type models. Although not well clari- 
fied, it has been strongly suggested by various experimen- 
tal measurements such as phase sensitive,^" tunneling^ 
ARPES^ penetration depthjiiili^ etc. that the pair- 
ing symmetry is also d wave for electron-doped cuprates 
at low and optimal doping. Thus one would expect that 
the rf-wave superconductivity in both electron- and hole- 
doped cuprates may be understood by a unified t-J type 
model. 

In this paper we study the t-t'-t"-J model by use of 
the slave-boson MF approach in the electron-doped case. 
In consideration of the AF order, the energy bands and 
the corresponding FS are calculated in the magnetic Bril- 
louin zone (MBZ). The experimentally observed FS evo- 
lution with doping is reproduced, essentially, a pocket 
appears first around (tt, 0), and another one shows up 
around (f, f) upon increasing doping. 

The t-t'-t"-J model Hamiltonian is written as 
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H = -t ^ (4c,, + h.c.) - t' (4c,. + h.c.) - i" J2 (4c,. + h.c.) 

{ij)a {ij)2<y {ij)3<y 

(ij) ia 



(1) 



where (), ()2, Os represent the nearest neighbor (n.n.), 
second n.n., and third n.n. sites, respectively, and the 
rest of the notation is standard. No double occupancy 
is allowed in the model. The Hamiltonian is essentially 
the same as that in the hole-doped case except here for 
electron-doping, t < 0, t' > and t" < OMiL Through- 
out the work |i| is taken as the energy unit. Since we 
mainly focus on the qualitative results, no model pa- 
rameters are fine-tuned. Typical values are adopted: 
t' = 0.3, t" = -0.2 and J = 0.3. 

We make use of the slave-boson transformation, i.e., 
CiCT = ^\fi<y with bi. bosonic holon operator, 
fermionic spinon operator and the constraint b\bi + 

J^rrfLfio- = 1 at each site. Since we are interested in 
low temperatures, boson condensation is assumed, i.e., 
(bi) = = y/S {5: doping concentration). Then we de- 
couple the Hamiltonian by defining the MF parameters: 

[thus m = (— l)*(5'f ) representing the AF order] and the 
uniform bond order 



X 



Based on the above treatment, the Hamiltonian can 
be expanded as follows in momentum space (up to irrel- 
evant constants) 

H = ^{ekfl^fka + Sk+qfl+Q^fk+Qa) 
k,cr 

-2JmY,<y{fljk+Qa + h.c.) + 2NJ{x^ + m^) , 

(2) 

where Sk — (2|i|5 — Jx){coskx + cos ky) — 
it's cos kx cosky — 2t"6{cos2kx + cos2ky) — fi, Q = (tTjTt), 
and N is the total number of lattice sites. The chemical 
potential is renormalized: fi — fiQ — X + J(l — S) with A: 
the Lagrange multiplier. (The local constraint has been 
relaxed to the global one.) Note that the summation to 
k is over the MBZ: — tt < kx i ky < tt. 

By use of the following unitary transformations: 



fka 
fk+Qa- 



COS 0k cr sin 0k 
a sin 0k cos 0k 



OLka 
Pka 



(3) 



with 



the Hamiltonian (O can be diagonalized into 
H = Y.^£.ko.al,aka + ^k^lM + 2iVJ(x' + vn?) , 

ka 

(5) 

with the energy bands 



ika,p = (£fc + £fe-HQ)/2 T \/(£fc+Q - £fe)V4+ (2Jm) 



(6) 



The free energy is simply given by (fcs = 1) 



F = -2T X Xln(l + e"«'"'/'^) + 27VJ(x^ + m2) . 

t)—a.,j3 k 

(7) 

The MF parameters m and x are then decided by 
(fm)^ ~ ^ ^"^^ (^) ^ "^^"^ chemical potential /i is 
adjusted to yield the right fiUing: = N{1 - 5). For 

each given doping 5 and temperature T, the quantities 
m, X and /i are calculated sclf-consistently. 
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cos 2ek = {sk+Q -ek)/\J {ek+Q - ekf + 4(2 Jm)^ , (4) 
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FIG. 1: The MF parameters m and x (negative) as functions 
of (5 at T = 0.001|t|. The model parameters are taken as: 
t' = 0.3, t" = -0.2, J = 0.3 (in units of \t\). 

In Fig.^we give the results for m and x as functions of 
doping at T ~ 0. It is seen that the staggered magnetiza- 
tion decreases with increasing doping, and goes sharply 
to zero at around 5 — 0.145. The AF order is overesti- 
mated due to the MF treatment: for example, in the un- 
doped case {5 = 0), one has m = 0.5, which is larger than 
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the accurate value 0.3 for the two-dimensional Heisen- 
berg model. In spite of this, the obtained energy bands 
are still instructive to understand the experimental ob- 
servation as shown below. 

The energy dispersions are plotted in Fig. [3 and 
the corresponding Fermi surfaces in Fig. |S| Note that 
the beginning Hamiltonian with t < (and t' > 
0, t" < 0) is already the transformed version after 
particle-hole transformation, thus doping electron ac- 
tually means doping hole in our treatment. For the 
same reason, ^ka,i3 in Fig- [21 should be inverted, i.e., 
S.ka,i3 —£,ka.p if they are understood as the en- 
ergy bands for the original electrons. The Matsubara 
Green's function for /-operator is simply written down: 
G{k,iujn) = cos"^ 9k/ {iujn - ^ka) + siv? Ok/ {iujn - £,kp)- 
And the single particle (c-electron) spectral function is 
A{k,uj) = -{5/TT)lvcLG{k,iuJn)\iuj„->u,+m- The density 
plots for J Auj A{k , uj) / {1 + e~^/'^) over an energy inter- 
val around the Fermi level have been shown in Fig.^for a 
few dopings. We come to see the details in the following. 

At low doping {5 = 0.04), the Fermi level crosses only 
the a-band. The low-energy states are around (tt, 0) 
and equivalent points, in consistence with the numer- 
ical results based on the single/slightly electron-doped 
systems i'^1^'^ Correspondingly small FS pockets appear 
around (±7r,0) and (0,±7r), as shown by the upper- 
left panel of Fig. |21 (see also the left panel of Fig. 0}. 
The energy gap above the Fermi level at (f, f ) is about 
0.44|t| ~ 0.15 eV (if typical \t\ = 0.35 eV is taken), which 
is in the same energy scale ~ 0.2 eV as measured by 
ARPES.^'^ With increasing doping {6 — 0.1), the stag- 
gered magnetization is reduced. Then the two bands be- 
come close to each other and the energy minimum of the 
/3-band at (f,f) shifts towards the Fermi level. The 
FS is still contributed by the single a-band: the pock- 
ets around (±7r, 0) and (0, ±7r) expand. Although the 
/3-band has not touched the Fermi level, it may con- 
tribute enough strong spectral intensity around (§, f) 
if it further approaches the latter with increasing dop- 
ing. This is evidenced by the middle panel of Fig. 0] for 
5 = 0.12, where finite spectral intensity is clear around 
(|-, ^). Also there, it is observed that half of the square 
around (tt, 0) loses much of its intensity, consistent with 
the ARPES findings. When doping continues to increase, 
e.g., for 6 — 0.144, the two bands are strongly over- 
lapped and both crossed by the Fermi level. New FS 
pockets centered at (±-|, ±-|) contributed by the /?-band 
will emerge, as shown by the solid lines in the lower-left 
panel of Fig. |3| If the FS curves for a-band (the dashed 
lines) are moved outside of the MBZ by the wavevector 
Q, one will see totally three FS pieces in the first quad- 
rant, which, when looked together, are close to a large 
curve around the point (7r,7r). This is more explicitly 
seen from the spectral intensity as shown in the right 
panel of Fig. ^ From the above, we can see that the FS 
evolution with doping is in essential agreement with the 
ARPES measurementsiS 

Eventually, when 5 = 0.145, the AF order breaks down. 



The two bands touch, which will merge into a single band 
when plotted in the original Brillouin zone. This is sug- 
gested by the lower-right panel of Fig. where a single 
continuous FS curve will form in the first quadrant if the 
dashed lines are moved outside of the MBZ by Q. 

The energy bands at finite doping as shown in Fig. [21 
are, if inverted, qualitatively similar to those obtained by 
Kusko et al. based on the t-t'-t"-U model, see Fig. 2 in 
Ref. . In order to realize the appearance of the pocket 
around (f , f) at a relatively large doping, Kusko et al^ 
introduced the doping-dependent U — Ueff(S). Then the 
issue, whether the pocket around (§,f) really appears 
or not, is dependent on the phenomenological form of 
Ueff{5). But the latter is not unambiguously clarified. 
In the current t-t'-t"-J model no parameters are intro- 
duced. Actually in our approach, the hoppings have been 
effectively multiplied by a factor S. We notice that the 
pocket around (-f, f) shows up naturally upon increasing 
doping. 

We would also mention that, due to the intrinsic lim- 
itation of the MF treatment, some ARPES data are not 
explained, for example, the development of the spectral 
weight inside the gap. This needs the study beyond MF, 
as recently done by Kusunose and Riceii& 

So far we have studied the FS evolution with doping 
based on the t-t'-t"-J model. Very recently the similar 
model was adopted by Li et ali^ to calculate the spin dy- 
namics in the SC and normal states around the optimal 
doping, and consistent results with neutron scattering 
data on Nd2-2;Cea;Cu04 {x = 0.15)M were obtained. It 
shows that the t-J type models may be available in the 
study of electron-doped cuprates. On the other hand, 
more theoretical predictions based on them are needed. 
As an example, following Li et al. one may further calcu- 
late the spin dynamics in the AF state. As already shown 
above, within certain doping region the two energy bands 
are both crossed by the Fermi level. It is expected that 
the interband excitation will lead to characteristic spin 
susceptibility. Another interesting topic is to investigate 
the potential spinon pairing under the AF background 
and the possible coexistence of the AF and SC states. 
Note that the pairing here is special due to the existence 
of two bands around the Fermi level. (Similar situation 
seems not to be present in the hole-doped case.) The 
studies on these issues are in progress. 

In conclusion, the electron-doped t-t'-t"-J model has 
been studied in the AF state. By using slave-boson MF 
treatment we have calculated the two energy bands in 
the MBZ. It was shown that at low doping only one 
band is crossed by the Fermi level and small FS pock- 
ets appear around (±7r, 0) and (0,±7r). Upon increasing 
doping, the other band will be crossed and new FS pock- 
ets form around (±-|, ±^). The evolution of the FS with 
doping and the calculated spectral weight are in good 
agreement with the ARPES data. 
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FIG. 3: The Fermi surfaces plotted in the MBZ: one-to-one correspondence to the energy bands shown in Fig.|5| 
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FIG. 4: The maps for spectral weight, obtained by integrating A(k,uj) times the Fermi function over an energy interval [—0.06, 
0.12]|t| (~ [—20, 40] meV corresponding to that adopted in the ARPES experiments) around the Fermi level. Highs are denoted 
by white and lows by black. 



